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Abstract 



We consider the asymptotic behavior of an evolving weakly coupled Fokker-Planck system of 
P I ■ two equations set in a periodic environment. The magnitudes of the diffusion and the coupling 

arc respectively proportional and inversely proportional to the size of the period. We prove that, 
as the period tends to zero, the solutions of the system either propagate (concentrate) with a 
fixed constant velocity (determined by the data) or do not move at all. The system arises in the 
modeling of motor proteins which can take two different states. Our result implies that, in the limit, 
the molecules either move along a filament with a fixed direction and constant speed or remain 
immobile. 

Key- Words: Hamilton-Jacobi equations, homogenization, molecular motor, singular perturbation, 
r**"" ! viscosity solutions 

. AMS Class. No: 35B25, 35B27, 49L25, 92C05 

^ ■ 1 Introduction 

We study the asymptotics, as e — > 0, of the weakly coupled Fokker-Planck system 



>< I ( n\ tt - sAn] - div^A^f )) + ^(i K = Ml H 



, 4 

S ■ { in M d x (0,oo), (1) 



Systems like ([T|) have been used to model the motion of motor proteins along molecular filaments 
or microtubules \12>\ [20} O [2Tj [91 [T3] • The intracellular transport in eukaria is attributed to motor 
proteins that transform the chemical energy into mechanical motion. For example myosins, which are 
known for their role in the muscle contraction, move along actin filaments and kynesins move along 
microtubules 

In ([T]) the molecules have the conformations 1 and 2 with densities n\ and n\ respectively and are 
influenced by the periodic potential ip provided by the filaments (see [H]). The (periodic) functions 
v 1 and v 2 indicate the rates of change between the two states. The existence of traveling waves, the 
asymptotic speed of propagation for large time as well as the presence of concentration effects for 
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models with diffusion and a periodic drift have been studied in many papers; for instance see [TOj [8]. 
A system similar to (pQ) also appears in the stochastic Stokes' drift where particles are suspended in 
a liquid and are subjected to diffusion and a net drift due to the presence of a wave in the liquid [7]. 
From the mathematical point of view, we refer to [I] for a related homogenization problem as well as 
to |15] for the homogenization of a similar equation with parabolic scaling. 

To formulate our result for the densities n\ and n 2 next we introduce the assumptions we will be 
using throughout the paper. In addition to 

v 1 ,v 2 and ijj are smooth and 1-periodic and u 1 > 0, u 2 > 0, (2) 

we assume that 



nl(x,0) > and n 2 £ (x,0) > for all x G R d , 

I® = nl(x,0)dx + / n 2 (x,0)dx -)■ I > 0, as e -)• 0, 
Jr Jr 

lim e _>o elnra^-jO) = lim £ _>.o £lnn 2 (-,0) = — oo locally uniformly in W~\ {0}, 
and 

limsup elnn* (y, 0) < for i = 1,2, 

E->0 

and, there exist constants A > and B such that, for all x G M. d and i = 1, 2, 

n*(x,0) < e^^-^l+B). 
Let 5 denote the Dirac mass at the origin. Our first result is 



(3) 
(4) 

(5) 



(6) 



Theorem 1.1 Assume 
sense of measures, 



), (0D, © and ©. There exists v G R d such that, as £ — > and in the 



n~.it, x) + n 2 (t, x) — v J(x — tv)Io- 



(7) 



To prove Theorem ll.il we analyze the behavior, as e — > 0, of the functions R\,R 2 : W~ x [0, oo) — > M 
which are obtained from n\ and n 2 by the classical exponential change of variable (we show later in 
the paper that n\ > and n 2 > 0) 



n\ = exp(— R e /e) and n 2 = exp(— R 2 /e), 



(8) 



and solve the system 



' R\ t - eA x Rl + \D X R\\ 2 - Z^(§ ) • + A^(f ) + ^ 2 (f) exp(^) = ^(f ) 

in M d x (0, oo) 

^ i? 2 1 - eA.U? + | JVfl 2 + ^(f ) exp(^) = v 2 (j ), 

(9) 

where we use the subscript x and ?/ to differentiate between differentiation with respect to x and 
the fast variable y = x/e. 

To state the second main result we recall that the "half-relaxed" upper and lower limits, denoted 
R l and R l respectively, of the family (R l £ ) e> o are given, for (x, t) G M d x [0, oo) and i = 1, 2, by 
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R l {x,t) = liminf R\(y,s) and Ft\x,t) = limsup Rl(y, s). (10) 

(y,s)-t(x,t) " (y,s)-+(x,t) 
e->0 e ^ 

We have: 

Theorem 1.2 Assume © and ©. T/ien 
(%) R 1 = Tf and R l = R 2 in R d x (0,oo). 

(ii) There exists a strictly convex H 6 C 1 (M ci ) satisfying H(0) = and, for some C > 0, H(p) > 
\p\ 2 — C such that R = R = R and R = R 1= R 2 are respectively subsolution and supersolution of 

R t + H(R X ) < and R t + H{FQ > mR d x(0,oo). (11) 

(mj Assume, in addition, (|5|). T/ien 

i? 1 = R 2 = R l = R 2 = +oo, in (R d \ {0}) x {0}, i? 1 (0,0)>0 and R 2 (0, 0) > 0. (12) 

('zwj The special direction v in Theorem lj.il is 

v = DH(0). (13) 

The asymptotic behavior of a time-independent version of (pQ) set in [0, 1], which is also controlled by 
the same effective Hamiltonian H, was studied in [27] where it was proved that, for asymmetric poten- 
tials, the mass concentrates at either x = or x = 1. The asymmetry condition of [27] is DH(0) ^ 0. 
This behavior is in agreement with our study of the time dependent problem. Indeed we prove here 
that, if v = DH(0) ^ 0, the proteins (mass) move (spread) with constant velocity v. Hence the mass 
concentrates, for large times, on one end point of the filament, if the latter is assumed to have finite 
length. We also refer to [25] for a large deviation approach for the asymptotic behavior of the sta- 
tionary solution of a similar model but with two potentials and to [26] for the study of flashing ratchets. 

By slight modifications of the proofs, all the results in this paper extend to systems with two 
potentials like 

n\ t - eAnl - div^D^f )) + ^(fK = Ml >l 

in R d x (0,oo), (14) 

nl it - eAn 2 £ - div x (n^^ 2 (f )) + ^ 2 (f )n 2 = \v\^ )n\, 

the only difference being in the value of the effective Hamiltonian H(-) -see Section[2j We also remark 
that we cal also consider without any difficulty systems with more than two equations. 

Our work is inspired from the ideas in wavefront propagation and large deviations [17} [5] , the method 
of perturbed test functions in homogenization [16] and the methods used in the study of the concen- 
tration effects [61 [M] and motor effects [2511271 [26] . 

Throughout the paper solutions are taken to be either classical, if smooth or, otherwise, in the 
viscosity sense. We refer to E] for a general introduction to the theory of the latter. In addition 
we denote by C positive constants which are independent of e but may change from line to line. 
Moreover, B r (x) is the open ball in M. d centered at x and of radius r > and B r (x) stands for its 
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closure. When x = we simply write B r and B r respectively. Finally many statements in the paper 
hold for i = 1,2 without any changes. Hence, unless necessary, we will not be repeating the "for 

1 = 1,2". 

The paper is organized as follows. In Section [2] we introduce the cell problem corresponding to ([9]), 
we recall that it has a solution and we introduce H. In Section [3] we present some preliminary facts 
about the family (re*) e >o and study the properties of the families (R £ ) s> q that are needed to prove 
the convergence to r\r l . Theorem 11.21 is proved in Section HI Using the results on the asymptotic 
behavior of the family (R 1 £ ) e> q we prove Theorem 11.11 in Section \E\ The asymptotic behavior of the 
family (re*) e> o in a more general setting is analyzed in Section [6l In Section [7] we compare our results 
with the case of the parabolic scaling in [15] . Finally, for the convenience of the reader, we present in 
the Appendix a sketch of the proof of the solvability of the cell problem and the properties of H. 

2 The cell problem 

In view of the presence of the exponential terms in ([9]) it is natural to expect that the Rl's and R 2, s 
converge, if at all, as e — > to the same limit R. Following [27] we insert in © the formal expansion 

Rl(x,t) = R{x,t)+e<i) l {^) + 0\e 2 ), 
and, keeping only the terms multiplying e° and writing y for the fast variable |, we conclude that 

R t - Ayf 1 + \D y 4> 1 + D X R\ 2 - D y 4> ■ (D y( j) 1 + D X R) + A y 4> + u 2 exp(<^ - 2 ) = u 1 

in R d x (0, oo) 

R t - A y <f) 2 + \D y (j) 2 + D X R\ 2 + v x exp((/> 2 - 1 ) = v 2 . 

(15) 

The goal is then to come up with 4> l, s so that (]15p is independent of y. This leads to the problem 
to find, for each p £ M. d , a unique constant H (p) such that the system, which is usually called the cell 
problem, 

- Ayf 1 + \DyCf) 1 + P \ 2 ~ Dy^ ■ [D y^ + P )+ Ayl/j + ^ CX?^ - ^ ) = ^ + H fr) 

in R d , (16) 

-A y 4> 2 + \D y 4> 2 + p\ 2 + v l exp((/> 2 - (j) 1 ) = v 2 + H(p), 

admits an 1-periodic solution ((f) l ,4> 2 ) called the corrector. 
We have: 

Lemma 2.1 For each p £ M. d there exists a unique constant H(p) such that f)16j) has an 1-periodic 
solution ((f) 1 , (j) 2 ). Moreover, H € C 1 (M C '), H(0) = 0, H is strictly convex, there exists a constant C 
such that H{p) > \p\ 2 — C , and , hence, H(p) — > oo as \p\ — > oo. 

A proof of Lemma 12.11 for d = 1 was included in |27] . For the convenience of the reader we sketch 
in the Appendix the proof in W 1 . 

When considering the system (|14|) with two potentials, (|16p is replaced by 

' -Aytf + \DyCj) 1 +P\ 2 ~ Dylf) 1 ■ (Dy^ + P )+ Ay^ + I/ 2 6X P (^ - ft) = ^ + H (j)) 

in R d . 

^ -A y 4> 2 + \Dy(j) 2 + P \ 2 - Dy4> 2 ■ {Dy(j) 2 + P )+ Ay^ + ^^{^ -^) = ^ + H \p) , 
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3 Some preliminaries and the properties of R % £ 

We summarize in the next lemma some basic properties of the families (n l £ ) £> Q. They are the con- 
servation of mass, the strict positivity of the n*'s, and a global upper bound yielding that, as e — > 0, 
there is very little mass at infinity. 
We have: 

Lemma 3.1 (i) For all t > and e > 0, 

nl(x,t)dx+ n 2 £ (x,t)dx = J°. (17) 



(ii) Assume ([2]) and ([3]). Then 

< n\ in R d x [0,+oo). (18) 
(Hi) Assume, in addition, ([6]). There exists D > such that, for all (x,t) el^x [0,+oo). 

,-A\x\ + B + Dt^ , . 

n\{x,t) <exp( 1 1 ). (19) 

In particular, there is small mass at infinity, i.e., for all t >0, there exists M = M(t) > such that 



L 



n\{t,x)dx — >0. (20) 

\x\>M 

Proof. The conservation of the mass follows from adding the equations in ([T]) and integrating over 



The form of ([T]) and (J2]) and ((3]) allow us to use maximum principle-type arguments to obtain (fT8|) 
and (fT9|) . 
Indeed let 

Fi tS (n) = n t — eAn — d\v x {nD y '4:{ — )) H — v (-)n v ( — )n, 

and 

i ? 2, £ (^) = nt — eAn H — v 2 {—)n 1/ 1 (— )n. 

It is easy to verify that n £ = min(n^,n^) and N £ = max(n^,n^) satisfy, in M. d x (0, oo), respectively 

max(F li£ (n e ),F 2j£ (n £ )) > and mm{F 1>£ (N £ ), F 2 , e (N e )) < 0. (21) 



Since is clearly a solution of the first inequality in (|2ip . (|18|) follows from the strong maximum 
principle. 

To prove (|19p we observe that, for D sufficiently large, 

G £ (j;, t) = exp( ) 

is a viscosity supersolution of the second inequality in (|2~T|) . Since, in view of ([6]), we also have 

iV £ (x,0) < G £ (x,0), 
we conclude using again the comparison principle. 
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Finally (|20p follows from (|19p after an appropriate choice of the constant M. □ 

We turn now to the properties of the i?*'s which are presented in Theorem 13.21 below. The proof 
is rather long. The first part, which provides an one-sided Lipshitz-type continuity in time, is based 
on the classical Harnack inequality. The lower bound in part(ii) follows from part(i). The arguments 
leading to the upper bound (part(iii)) are more tedious and require as an intermediary step, namely, 
the construction, again using part(i), of an appropriate local upper bound. 

We have: 

Theorem 3.2 (i) Assume ([2]) and ([3]). For all 5 > 0, there exists Cg > 0, such that, for all e > 
0, \z — z'\ < e, e5 < to and i,j = 1, 2, 

Ri(z',t + e)-Ri(z,t )<eC 5 . (22) 

(ii) Assume ([3]), and (JSJ). For any a £ (0,oo), there exists Eq = £o(a) > such that, for all 

R £ > -a in R d x [0,oo). (23) 

(Hi) Assume ([2]), ([3]), (HJ) and ©. For any compact subset K ofM, d x (0, oo), there exist Ck > 
and £\ = E\{K) > such that, for all e < E\ and (x,t) £ K, 

Rt(x,t)<C K . (24) 

Proof. We begin with the 

Proof of (I22D : Observe that n l (y, r) =n^(ey,er) and n 2 (y,r) = v? £ (Ey, et ) are positive (recall (]18p ) 
solutions to 

n]. — AyU 1 — n 1 A y ip(y) — D y ip(y) • D y n l + v l n l = u 2 n 2 

in R d x (0,oo) (25) 

r? T — A y n 2 + v 2 n 2 = v l n l , 

a linear parabolic system with bounded, according to ([2]), coefficients. It follows from the classical 
Harnack inequality [23] that, for each 5 > 0, there exists, an independent of e, Cs > such that for 
all y e R d , t >5 and i, j = 1, 2, 

sup n i (z,T )<C s inf n j (z,T + l). (26) 

Rewriting (I26p in terms of n 1 and n 2 and in the original variables (x, t) we get 

sup n l £ (z,to) < Cs inf n £ (z,t +e), for (x ,to) e^x [ed,+oo). (27) 

«e5 E (i ) zeB £ (x ) 

Finally using ((8|) we obtain (|22|) . 
We continue with the 

Proof of the uniform bounds from below: Arguing by contradiction we assume that for some 
a > there exist e^ — > and (xf.,tf.) € R d x [0, oo) such that 

min (Rl k (x k ,t k ), R 2 k (x k ,t k )) < -a. 
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Newt observe that, in view of © and (J6]), for s < e a with e a small enough, we have 

a 
2' 



mm(Rl k (x k ,t k ),Re k (x k ,t k )) > 



Moreover min (i?*,i? 2 ) is a supersolution to 

max (Rl t - eA x Rl + \D X R\\* - £>^(§) • D,^ 1 + A^(f ) + ^ 2 (f ) - ^(|), 
JS 1 - eA.i? 2 + |A,U?| 2 + ^(f ) - ^ 2 (f )) > o, 



which admits — | — ct as a subsolution provided c is chosen sufficiently large. 
It follows that, for e < e a and t > 0, 



min (Rl(-,t),R 2 £ (-,t)) > ~ - ct, in R d , 
and therefore, if 5 = for e < e a , 

min {Rl,R 2 £ ) > -a, in M d x [0, 5}. 
As a result the sequence (x k ,t k ) chosen at the beginning of the proof must satisfy t k > 5. 



Using ([22]) we deduce that there exists Ci > such that, as k — > oo and for all cc such that 
\ x — x k\ < £k and i = 1,2, 

R\ k {x, t k + e fc ) < Ri k (x k ,t k ) + < -a + Ci£fe. 

It follows that 



f . f 2£fc!_!!^I!^: -2~ci 

/ n l £ Jx,t k + e k )dx > / e " e fc cte > |S £fe (x fe )|e e fc 



The right hand side of this inequality blows up as k — > oo, while the left hand side is bounded in 
view of (fTT]) and (j3|), again a contradiction. 

The last part of the proof is devoted to the 
Proof of the uniform upper bounds on compact: Fix a compact subset K of M. d x (0, oo), 

observe that 

to = inf{s G (0, oo) : there exists x € M. d such that (x, s) € K} > 0, 

choose t\ S (0,to) and write i\ = ti/2 and £2 = ii/4. 

It follows from (TTTj) and ([20]) that there exist e\ > and M > both dependent on t\ such 
that, if e < E\, then 

/ n*(:c,ii) + ra 2 0,ti) dx > 

./|x|<Af 1 

and, hence, in view of ([8]), there exists some a > such that 

min i^O, ti) < 6= -ln(^). 
|x|<M £V ' ~ K M d ' 

1=1,2 



7 



Assume next that the above minimum is achieved at some point x £ € Bm and for i = i* . Applying 
([22]) L = |_^J times with i = i* , j = 2, xq = x £ and S = i\, we find some C = Cs > such that, for 
all x G B Le (x £ ), 

R 2 £ (x, i 1 +Ls)<b + CLe <b + Ci X - (28) 



Choose 7 > b + Cii and for some /3 > to be fixed below define (ft 1 : -B^ (x e ) x (0, oo) — > R by 



1 



^ 1 ^ t ) = f2_| T _ T ,2 +^ + 7 

We claim that, for e < £2 = rnin(i2,£i), 

i? 2 < 1 in Q £ = B h (x e ) x [t! + Le, +00), (29) 

and, therefore, 

R £ < <ft l in Ql = B h {x e ) x [ti,+oo). (30) 
To prove (129p we first notice that, in view of (128|) and the choice of 7 

i^(-,t 1 +Le)<& + Cti<7<^ 1 (-,ti + Le) in B tl (x £ ) 



Moreover, if (3 is large enough, using ([2]), for all (x,t) £ Q e , we have 

a ,1 m ,12 « / 2d 8\x - x £ \ 2 \ 4\x - x £ \ 2 ^ 2 , . 

-A^ + |^| = P-e [ {n _ {x _ Xe{2)2 + ^ _ {x _ Xg|2) 3 J + (tf - |, - ge |2)4 > g = ma -^ ( 

The inequality above and ([9]) yield that, for all (x,t) £ Q £ , 



<ft\ - £A x (ft 1 + \D x (ft l \ 2 + exp( ^ - Rl£ ) > v 2 . 

Since R 2 < ft 1 on the parabolic boundary of Q £ , (|29p follows from the maximum principle. 
For e < £2 let Q 2 = (R \ Bf 2 (x £ )) x (ii,+oo) and for positive constants a, rj, £ to be fixed below 
consider the map <^> 2 : (R d \ Bf 2 (x £ )) x (ti, +00) — >■ R given by 



</> 2 (x,t) = ' /' +vt + C- 

t — t\ 

We claim that 

R 2 < (f> 2 in Q 2 . (31) 

As above we will show that, if we choose a, r\ and C appropriately, (ft 2 is a supersolution in Q 2 of the 
equation satisfied by R 2 and is above R 2 on the parabolic boundary of Q 2 . 
To this end notice that, in view of (|30l) . we may select £ and 77 large enough so that, for e < £2, 

R £ < r]t + ( < eft 2 on &B fa (x £ ) x (t x , 00). 

Moreover, for sufficiently large a, 77, we have 

,1 . ,1 , ^ ,19 a\x — x £ \ 2 2£ad Aa 2 \x — x £ \ 2 „ 

»,'-A^ + | D ^=,--l-^i-- — + ' ' >g, 



2 
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and, hence, in Q 2 , 

Since clearly R 2 {-,h) < ^{'M) in (M d \ B h {x £ )), using again the maximum principle we find (131 p . 

To conclude observe that (I30j) and (|31|) yield that the family (R 2 ) £> q is uniformly bounded from 
above in any compact subset of W 1 x (ti,oo[ and thus, in particular, on K, for e < £\. Finally, 
using again (I22p . we deduce that the family (R\) £> q is also uniformly bounded from above on K, for 
e < min(ei,t - h)- □ 

4 Convergence to the Homogenized equation-The proof of Theo- 
rem 11.21 

Before we begin the proof, we remark that, in view of the claimed properties of H (convexity and 
coercivity), the equation 

R t + H(R X ) = in R d x (0,oo), (32) 

admits a comparison principle even for initial data taking the value oo (see [32]). However, this and 
Theorem 11.21 do not lead to R < R. Indeed we show in the next section that, in addition to (|12p in 
Theorem 11.21 the conservation of mass yields 

E(0, 0) = ^(0, 0) = R 2 (0, 0) = 0. 

l 2 

To be able to use the comparison principle it is necessary to have R (0, 0) = R (0, 0) < 0, which is not 

1 2 1 2 

possible. Indeed, since R and R are upper semicontinuous, (fT2|) implies R (0, 0) = R (0, 0) = +oo. 
We continue with the 

Proof of Theorem II. 2L It is immediate from (I22p that (i) holds. Moreover it follows from Theorem 
EOl (ii) that 

< R 1 < R 1 and < R 2 < R 2 on R d x [0, oo), (33) 

and thus, in particular, the second part of (fT2|) . 

To prove (jlip we employ the so called perturbed test function method. Since the arguments are 
similar here we only show that R is a subsolution of (lll|) . 

To this end we assume that, for some smooth ip } R — <p achieves a strict local maximum at 
(xo,£o) £ ^ d x (0, oo). Following [27], we perturb ip using the solution ((f) 1 , (f) 2 ) of the cell prob- 
lem (fT6|) corresponding to p = Dtp(xo,to). It follows that there exists (x £ , t £ ) S M. d x (0, oo) — > (xq, to), 
as e — > 0, such that the max^i^-R^ — tp — ecj) l (-)) is attained at (x £ ,t £ ) and, without any loss of 
generality since the argument is identical, for i = 1. Hence, 

Rl(x £ ,t £ )-e^(^) > RKxM-e^ 2 ^). (34) 

That R\ is a solution of Q yields 

(ft(x e ,t E ) - £A x p(x £ ,t £ ) - A y ^(f) + \D x <p(x e ,t E ) + Dy^if )\ 2 
-D„V(f )• (DM^,t e ) + D x <f>\^ )) +A^(f ) + ^ 2 (f )exp(%^) <^(f ). 
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Using the latter and (|34p and writing p £ = Dip(x e ,t e ) we get 

<p t (x £ ,t £ ) - eA xV (x £ ,t £ ) ~ A^Hf) + \Pe + D y ^(f)\ 2 - A^(f) • (p c + A^(f )) 
+A,V(f ) + ^ 2 (f) exp(^(f ) - ^ 2 (f )) < ^(f )• 



It follows from the above inequality and the definition of the effective Hamiltonian (|16j) that, for 
some o(l) — > as e — > 0, 

<p t (x e ,t s ) - eA x (p(x e ,t E ) + H(D x tp(x £ ,t £ )) < o(l), 

and the conclusion follows letting e — > 0. 

To prove the first part of (|12p . since the arguments are identical, here we only show that R}(-, 0) = 
R (-,0) = 00 in M rf \{0}. To this aim, we fix k > 0, we select an auxiliary function £ £ C°°(IR d ) such 
that < £0) < 1 for i£K d \ {0} and £(0) = 0, and prove that 

max(Rj + H(lg.),I£ - k£ ) > in M d x [0, 00). (35) 

Since we already know that (I35p holds in R d x (0,oo), to conclude we assume that, for a smooth 0, 
i? 1 — achieves a (strict) local maximum in (xq,Q) and we prove that either 

^(zo.O) > k£(x ), 



or 



<Mx ,0) + tf((fe(xo,0)) >0. (36) 
If xq = 0, the former is clearly true since ^(0,0) > = «f(0). 

So we assume that i? 1 — has a local maximum in (xo,0) with xq 7^ and, in addition, that 
R (xq,0) < k£(xq). Repeating the arguments used earlier in the proof we obtain, for some (x £ ,t £ ) — >• 
(xq, 0) as e — > 0, we have 

(p t (x £ ,t £ ) - eA x (p(x £ ,t £ ) + H(D x ip(x £ ,t £ )) > o(l). 

Indeed §5§ and the facts that ^(xo, 0) < k^(xq) < k and lim^o R\{U: 0) = +°° f° r an V near yield 
that t £ > 0. The claim now follows by letting e — > 0. 
Assume next that, for some xo 7^ 0, ij^xo, 0) = b < +00. We fix 5 > and let 

I |2 
(5/ .\ X — Xo 

A* 0=,*) = j 7* 

for 7 = 7(5) > to be determined later. 

Since R} is lower semicontinuous, i? 1 — /i 5 attains a minimum at some (x^,^) Gl^x [0,oo) such 
that, as S — > 0, (xs,t$) — > (xo,0), provided that 7 — > 00 as 5 — > 0. 

Observe next that the choice of (x^,^) yields 

1X5 ~ X ° |2 < fl 1 ^, *<s) + |X '~ X012 + 7*«5 < R\x Q , 0) = 6, (37) 
d 

and, hence, 

|x5 — xo| < Vbd. 
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If ts > 0, according to part (ii), we must have fj,^(xs,tg) + H(Dfj,°(xg,ts)) > and, hence, 

-1 + H(- 2{XS ; X0) )>0, (38) 



which cannot be true if we choose 7 > sup^ K ^ H(—^-). 

Now we assume that ts = 0. If R}(xo,0) < ^(^0)) then ([3T|) yields R}(xs,0) < ks£(xs) for some 
k$ — > k as 5 — > 0. Using ([35]) with n$ at the point (xg, 0) we obtain again (|38|) and thus a contradiction. 

It follows that ^(xojO) > k£(xo), which also leads to a contradiction, since it holds for arbitrarily 
large k and R}{x§, 0) = b < +00. 

The first part of ()12|) now follows. □ 



5 The transport of the concentration points - The proof of Theo- 
rem 11.11 

We present here the 

Proof. [Proof of Theorem 1 1.1] Using (jlip . (|12p . the standard optimal control formula [22 ] 1181 [3] and 
a barrier argument similar to the one of Section 5 in [17], we obtain that R satisfies 

R(x,t)> inf / F*(C(s))ds + max(i? 1 (0,0), ; R 2 (0,0)), (39) 

(C( s ), s )eR d x(o,«D)+ Jo 

<(0)=0, C(t)=x 

with #*(p) = sup 96Rd (p • q - H(q)). 

Observe next that, since Lemma 2.2 yields 

lim p-q — H(q)> lim p ■ q — \q\ 2 + C = — 00, 

\q\— >+oo \q\— >+oo 

the maximum of (p ■ q — H(q)) is attained at some q p G M d such that p = DH(q p ). 
Moreover, since H is strictly convex and H(0) = 0, we have 

DH(q p ) • q p — H(q p ) > with equality only if q p = 0. 

Hence we deduce that 

H*(p) > for all p / DH(0) and H* (DH(0)) = 0. 

Therefore, using (]39h and (]12H . we obtain that 

# > in M d x M+ and if R(x, t) = 0, then x = tDH(0) and ^(0, 0) = R 2 {0, 0) = 0. (40) 

It follows from the latter, © and (fT0|) that, for i = 1, 2, 

limsup <(y,s) = in R d x [0, +00) \ {(i£>lf(0), t) 1 1 G [0,+oo)}. 

{y,s)^{x,t) 

e->0 

Finally the last claim and (fT7|) yield that the {n\ + ra^'s converge weakly along subsequences to a 
measure n with 

suppnC {(tVl?(0),i) |t G [0,oo)}. 
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Since we also know that, according to (f20|) . no mass escapes to infinity as e — > 0, we deduce (|7|) 
using (I17j) . Moreover, in view of (140p . we obtain that 



^(0,0) =i? 2 (0,0) = 0. 

□ 

6 The case with several Dirac masses initially 

We proved Theorem 11.11 under assumptions © which imply that the densities n\ and n 2 . are both 
initially concentrated at the origin. The result can be generalized to densities concentrated at several 
points and probably not on the same points. If this is the case, the initial condition is written as 

lim e _5.o elnn^(cot,0) = lim e _s.o £hin^(-,0) = — °o locally uniformly in R d \ (AuB), 
and / 41 j 

lim sup e In n\ (y, 0) < for all x 6 A U B. 

y—tx 

K £->0 

We have: 

Theorem 6.1 Assume ([2]), (|3j), (jH), © and (j4T|) . Then, for i = 1, 2 and as e — >• 0, along subse- 
quences and in the sense of measures, n\ — 1 n 1 and n 2 . — n 2 with with 

supp (n 1 + n 2 )(-, t) C C(t) = {xi + tDH(0), ■ ■ ■ , x n + tDH(0)} U {y x + tDH(0), ■ ■ ■ , y m + £1)17(0)}. 

Proof. The proof of Theorem 16.11 follows along the same lines as the one of Theorem 11.11 The only 
difference is that (fT2|) and ([39)) are replaced respectively by 

R 1 = R 2 = R 1 = R 2 = +oo in (R d \ C(0)) x {0}, R 1 > and R 2 > in C(0) x {0}, 

and 

R(x,t)> inf j / ^(C(s))ds + max(i? 1 (C(0),0),/i 2 (C(0),0)) 1 . (42) 

(C(s),s)m d x(0,oo) [Jo J 
C(0)eC(0),C(t)=z 

Therefore, we have 

1J>0 in M d x(0,oo) and if lj(x, t) = 0, then xeC(t). 
The other parts of the proof are similar. □ 

7 A comparison with results for parabolic scaling 

We describe here the connection between our result and the study in [15] of the asymptotics, as e — > 0, 
of the solutions to 

n e ,t = (a>ij(-, -n)n E , Xj )xi + -n)n e , Xi + -»c(-, -^)n £ in R d x(0,oo). (43) 



12 



Notice that this is a single equation — not a system — obtained after a parabolic scaling (x, t) h-» 
(f,4r) rather than the hyperbolic scaling i-)- (|,|) of the problem we consider in this paper. 

Nevertheless, as we explain below, there are some similarities. 

It is proved in [15] that the solution n £ of (|43p admits the expansion 

n e (x,t) = w(-,\)exp(-^)v°(x--t,t) + o(l). (44) 
The result we obtain here with the hyperbolic scaling formally gives 

nl(x, t) = w\^, (i)<5(z - M(O)t) + o(l). (45) 

In particular we obtain a Dirac mass instead of the function v° while the term exp(— ^) disappears 
because we have conservation of mass. Moreover t/e is replaced by t, because of the difference in the 
scaling. There is, however, as we explain below a close connection between the function w in (|44j) and 
w i in (fi5|) . 

Indeed iu in [15] is the principal eigenvector of the periodic cell problem 



w s - (aij (z, s)w Zj ) z .-bi(z,s)w Zi - c(z, s)w = A w, 

while we have 

w % = exp(-</> 1 ), 

where (cp 1 ,^ 2 ) is the principle eigenvector for the cell problem (|16p corresponding to p = with the 
corresponding eigenvalue H(0) = 0. This is because for (x, t) in the support of n*, we have R(x, t) = 
and thus DR(x, t) = 0. 

A The proof of Lemma 12.11 

Step 1: We prove that, for all p S M. d , there exists a unique H(p) with the properties stated in Lemma 
EH Following [27], for i = 1, 2, we define 



(46) 



X l (y) =exp(-p-y-cj) i (y)). 
Multiplying the two equations in ()16|) by X 1 and X 2 respectively, we obtain 

-AyX 1 - div y (DyiJ; X 1 ) + v x X x - v 2 X 2 = -H(p)X 1 , 
-A y X 2 + v 2 X 2 - v x X x = -~H(p)X 2 . 

with the boundary condition 

y -> e p - y X l (y) is one-periodic and X* > 0. (47) 
We also impose the normalization 

\x\y) + X 2 {y))dy = l. (48) 



o 
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It follows from the Krein-Rutman Theorem that, for all p € K d , there exists a unique constant H (p) 
and functions (X^X 2 ) satisfying (J46]) together with (|47|) and (f48l) . 

Step 2: To prove that -ff G C 1 (M d ) we rewrite ([16]) in terms of w p (y) = e p ' y X l (y), which in view of 
(|47p are 1-periodic, and obtain the new system 



-Aywl + 2p ■ Dywl - div y (Dyifj w^) + (~\p\ 2 +p ■ D y ^ + u l ) - v 2 w 2 p = -H(p)wl, 

(49) 

-AyWp + 2p ■ D y w 2 + {-\p\ 2 + z^ 2 ) u; 2 - v l wl = -H(p)wl. 

Assuming for the moment that Wp and w 2 are differentiable with respect to p, after differentiating 
491) we get 

-AydpWp- + 2p • Dydpwl - div y (Dyip dpW l p ) + (-|p| 2 + p • Vyip + v l ) dpw l p - v 2 d v w 2 
+2D y wl + (-2p + D y ip)wl = -H(p)d p wl - H'(p)wl, 



-AydpW 2 + 2p • DydpW 2 + (-|p| 2 + z^ 2 ) dpW 2 - v 1 d p Wp + 2A,u> 2 - 2pw 2 = -H(p)d p w, 
-H\p)w 2 p. 



Let {Wp x, w 2 ^) be the solution to the adjoint system of (jMD -note that Fredholm's alternative implies 
the existence of such a solution. Multiplying the equations of (j4T)j) by u;* and u> 2 *, integrating with 
respect to y and adding the two resulting equations we obtain 

2 J w l,* D v w ldy + 2 J w 2 ^D y w 2 dy + j DyipWpW^dy = (2p - DH(p))(j wlw^dy + J w 2 pW 2 pif dy), 

which yields a formula for DH(p). 

To prove the above claim rigorously, we write the system satisfied by difference quotients with respect 
to p and we use the same idea as above to prove that the difference quotients h" 1 (H(p + he^) — H(p)) 
converge to the above formula, as h — > 0. Indeed from this formulation we first obtain that H(p) is 
continuous with respect to p. Then we show that Wp and w 2 are continuous with respect to p and, 
finally, we pass to the limit h — > 0. 

To prove that Wp and w 2 are continuous with respect to p, we show that, for any sequence (p n )neN 
such that p n — > p, the corresponding eigenfunctions (Wp n , w 2 n ) converge to (Wp,w 2 ) as n — > oo. This 
follows, for example, from the stability of viscosity solutions, the continuity of H, the uniqueness (up 
to a multiplicative constant) of the eigenfunctions (Wp,w 2 ) and the normalization condition (|48p. 

Step 3: The proof of H(0) = 0. 
The adjoint system to (}4"6j) is 

-A v u l + D v ib • D v u l + v l u l = v l u 2 - UMu 1 , 

_ (50) 
-A y u 2 + z^ 2 u 2 = z/ 2 ^ 1 — H(p)u 2 , 

with the condition 

y — > e~ p ' v u l {y) 1 -periodic and u % > 0. 
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When p = 0, (|50|) admits the trivial solution (tt 1 , u 2 ) = (1, 1) and the constant H (0) = 0. The claim 
follows from the uniqueness of the positive eigenvector. 

Step 4: The proof of the strict convexity of H(p). 

Arguing by contradiction we assume that there exist pi,P2 such that 

Let (4>\,(p1) and (02j0|) be solutions of (fl~6j) corresponding to p = pi and p = p 2 respectively. 
Adding the corresponding equations we find 

+ l/ 2 exp( ^iz^l) 

< -Ay(^k) + \\D y cp\ + Pl \ 2 + \\D y 4>\ +p 2 ? - D y iP • (^4±El) _ ^ . ( £^±vi 
+A y iP + *f exp(</>{ - 0f) + % exp($ - 2 ) 



and 



-A„(&±&) + |A,(^) + £^| 2 + ^ 1 exp( ^ + ^"^-^ ) < v 2 + H(E^) 

' 4>\+<t>l <t>l+4>l ■ 



It follows that the pair ( * 2 , 1 2 2 ) is a strict subsolution to the cell problem f)16|) corresponding 
to p = 2i+£2_ This, however, contradicts the fact that H( Pl ^ P2 ) is the principal eigenvalue of the 
system. 

Step 5: The proof of H(p) > \p\ 2 - C. 
Rewriting (fl~6j) we find 

f -Aycj) 1 + |ZV>T - (D„^ - 2p) • Dycj) 1 -D y ^-p + Ayij + z^ 2 exp(^) 1 - </» 2 ) = i/ 1 + - \p\ 2 , 

\-Ay(j) 2 + \Dy(j) 2 \ 2 + 2 P ■ Dy^ + ^ 1 6X P (0 2 ~ ^ ) = if +H(p) ~\ P \ 2 . 

Assume next that maxr 0i i](^> 1 , 4> 2 ) is attained at the point y. If 4> l {y) > </> 2 (y), then 

-A(f> l {y) > and D^{y) = 0, 

and, hence, 

H{p) > -\\Ai/>\\ L ac - \p\ \\D y ijj\\ L oo - H^Hl- + \p\ 2 - 
If 4> l {y) > 4> 2 {y), then, similarly, we find 

h(p) > -II^IIl- + \p\ 2 . 

The claim now follows. 
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